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Abstract 



Passive scalar dynamics in wall-bounded turbulence is studied via Di- 
rect Numerical Simulations of plane channel flow, for a friction Reynolds 
number i?e» = 160 and a Schmidt number Se — 1. Peculiar to the present 
research is that the spatial resolution reaches far beyond what has been 
^ ■ employed in similar past studies. Our aim is to examine the statistics of 

f| I ' the most dissipative events across the various layers of the channel flow, 

jyj . and to compar e them to the homogeneo us is otropic case, where t he re - 

Q ■ cent studies by ISchumacher et al\ (|2005l ) and IWatanabe fc Gotohl foOOJl ) 

' TT* i have described a range of scalar micro-scales that require extremely high 

^►y spatial resolution to be properly resolved. 

^ r*| Resolution effects are observed on integral-scale quantities such as the 

O | mean profiles of the scalar dissipation and its variance. By examining 

probability distributions, it is found that the finest resolution is essen- 
tial for correctly computing small-scale statistics in the near-wall region 
of the channel. As expected, this high-resolution requirement extends 
w-v ' outwards to the channel centerline, where the behavior of isotropic tur- 

f — | bulence is recovered. However, high-intensity scalar dissipation events 

f~^. . are overemphasized by marginal resolution near the wall, while they are 

C*") ' underemphasized in the central region. 

in 

q ; 1 Introduction 

A passive scalar in a fluid flow is a diffusive physical quantity whose dynami- 
cal effects on the velocity field can be neglected. Examples of quantities that 
j^ ■ can often be conveniently treated as passive scalars are pollutants in liquids 

J~j | and gases, moisture in air, temperature in heated fluids (provided temperature 

fluctuations are sufficiently small), dye in liquids, chemical reactant species in 
industrial processes, etc. A comp rehensive review of the literature on passive 
scalars and turbulence is given by I Warhafq (J2000t . 



The turbulent mixing of a passive scalar arises from two different yet concur- 
ring mechanisms: advective dispersion by large-sca le motions, and small-scale 
diffusive molecular mixing. E xperimental results (I Warhaftl . 120001) as well as 



Direct Numerical Simulations ( Brethouwer et all 120031 ) show that these two 



physical processes are strongly coupled. As a consequence, a precise knowledge 



of the scalar's smallest scales is required for studying its integral-scale prop- 
erties, such as the mean profile, the transport, the production and dissipation 
rates of scalar variance, and the scalar flux. 

Understanding the smallest scales of scalar motion plays a key role in mod- 
eling the budget equation for the scalar variance, which is of interest in the LES 
and RANS approaches to turbulent flows. A satisfactory knowledge of the small- 
est scales mi ght also improve our und erstanding of turbulence as a whole, as 
discussed by IShraiman fc Siggial (|2000t ). Despite the linearity of the advection- 



diffusion equation with respect to the scalar, considerable internal i ntermittenc^ 



and s mall-scale anisotropy charact erize the passive scalar field (jSreenivasa: 



icy 
in, 



199 It ISreenivasan fc Antonial . 119971 ) . and are lar g er th a n the ones of the cor 



responding turbulent velocity fiel d ( Chen fc Cao , 19971: Mvdlarski fc Warhaftl . 
1998HWang et all . I1999J: fWarhaftl . l2000t lYeung et aZj . l2005h . 

An effective tool for the numerical study of the smallest features of scalar 
motions is the Direct Numerical Simulation (DNS) of the Navier-Stokes equa- 
tions. To obtain reliable and accurate results, however, an adequate spatial 
resolution must be chosen in order to resolve all the significant scales of motion. 

The typical lengthscale of the smallest scalar motions is the mean Batchelor 
length scale (tje), given by the mean Kolmogorov length scale (77) divided by 
the square root of the Schmidt number Sc — v /■y: 
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In the expressions above, v is the kinematic viscosity of the fluid and 7 is the 
scalar diffusivity; the operator (•) denotes an average taken over homogeneous 
directions and time. However, the fact that, on average, the smallest scales of 
the scalar motions are given by (77s) does not imply that a spatial resolution 
set at (tjb) is sufficient. 

The prope r spatial resolution f or the DNS of passive scalars has been r ecently 
addressed by ISchumacher et all (|2005| ) and IWatanabe fc Gotohl (J2007I ) in the 
context of homogeneous and isotropic turbulence. The general conclusions of the 
two papers are similar: they point to the existence and dynamical relevance of 
scalar micro-scales that are misrepresented, should one set the spatial resolution 
of the numerical simulation based on (tjb)- In fact, they demonstrate how the 
value of 7]b may become very small locally, and establish that in such cases scales 
larger than t\b and smaller than {tjb) must be represented in the numerical 
simulation for th e small-scale statis t ics to be correctly predicted. 

In particular. ISchumacher et all ( 2005T ) examined how these micro-scales in 
homogeneous turbu l ence b ehave as the Schmidt number is increased, whereas 
Watanabe fc Gotohl (|2007l ) used the computing power of the Earth Simulator 
computer to carry out the same analysis at increasing Reynolds numbers with 
fixed Sc = 1. The present paper deals with the same issue, and our original 
contribution will be to extend the analysis to wall-bounded turbulent flows. 
The plane channel flow is considered as the simplest prototypical wall flow that 
offers at the same time a simple geometry and the physical features of more 



complex wall-bounded flows. Aim of the paper is to show that the spatial 
resolution adopted in past DNS studies of wall turbulence with passive scalar 
does not resolve an important range of very small yet dynamically significant 
micro-scales, which are responsible for the extremely intermittent nature of the 
scalar field. 

To demonstrate this statement, three DNS of the same channel flow at differ- 
ent spatial (and consequently temporal) resolutions are carried out. The lowest 
resolution corresponds roughly to the standard spatial resolution adopted in the 
past for most passive scalar studies: it is of the order of a few wall units, and 
comes from using for the passive scalar the spatial resolution typically employed 
for the velocity field. This spatial resolution is larger than (ijb)- Our most re- 
solved simulation, on the other hand, possesses a spatial resolution that goes 
well below the (tjb) level. The effects of varying resolution on the statistical 
features of the scalar field can thus be addressed, and their change with the wall 
distance described. 

In particular, the focus of the present paper will be on the correct repre- 
sentation of the scalar dissipation. The dissipation eg of the scalar variance is 
defined as: 

A fde'V 

%— i x * 

where 8' = 8 — (8) is the fluctuation of the passive scalar about the local mean. 
Of course eg is particularly sensitive to the s mallest scales of motion: it is known 



that eg possesses very fine spatial structure (jSchumacher et all 120051), high in- 
ternal intermittency, and significant small-scale anisotropy ( Warhaftl . 120001 ). In 
addition to being a suitable indicator of the smallest scalar motions, eg is sig- 
nificant in many physical processes which are central in both industrial and 
environmental fields. A correct modeling of both integral- and micro-scale 
scalar dissipation features can, for example, establish whether a chemical re- 
action occurs or not, and allows us to correctly predict quantitative features of 
environmental pollution aimed at establishing appropriate safety thresholds. 

The outline of the paper is as follows. In fj2]the characteristics of the numeri- 
cal method will be briefly outlined, and its validation will b e discussed in £|2.1lbv 



reproducing the results of a channel flow DNS performed bv lJohansson fc Wikstrom 
(|1999h . In ^3] the numerical parameters employed in the simulations will be in- 
troduced, with a view to comparing between the spatial resolutions used here 
and the typical resolution employed in the past. The main results will then be 
presented, illustrating first in £0] the effects of the spatial resolution on integral- 
scale scalar quantities (mean profile and variance of scalar dissipation) and then 
in fj5] the effects of varying resolution on the micro-structural features of the 
scalar field. Finally, a brief summary and some conclusive remarks will be given 

in ga 



2 The numerical method 

The DNS code employed in this work has been developed f rom the pseudo- 
spectr al, mixed-discretization, parallel algorithm introduced bv lLuchini fc Quadrid 



(|2006r ) for the DNS of the velocity field for wall-bounded turbulent flows. The ex- 
tension of the original code to include the dynamics of a passive scalar has been 
implemented with the key requirement of keeping the same parallel strategy, 
and thus retaining the same computational efficiency. Since the wall- normal ve- 
locity - wall-normal vorticity formulation is used for the momentum equations, 
the evolutive equation for the passive scalar: 

m_ m__ a 2 e 

dt % dxi dxidxi 

is written to be formally identical to the wall-normal vorticity equation. 

We indicate with x, y and z the streamwise, wall-normal and spanwise direc- 
tions respectively. The corresponding velocity components are u, v and w, and 
the passive scalar is 9. The computational domain has extensions L Xl L y = 26 
and L z in the corresponding directions. The friction Reynolds number Re* is 
based on the channel half-width S, on the fluid's kinematic viscosity v and on 
the friction velocity u*, and is defined as Re* = u*8/v. The scalar concentration 
is expressed through the Schmidt number Sc = is/"f. 

A Fourier expansion is used in the homogeneous directions, whereas fourth- 
order accurate, explicit compact finite-difference schemes are used to compute 
derivatives in the y direction. The number of discretization modes (points) is 
indicated with N x , N y and N z . The collocation points in the near- wall direction 
are smoothly stretched from the wall to the centerline. Full dealiasing is used in 
the homogeneous directions. The time-integration algorithm uses a third-order 
Runge-Kutta scheme for the computation of the convective non-linear terms and 
a second-order Crank- Nicholson scheme for the evaluation of the viscous-implicit 
terms. 

The flow in the channel is made statistically stationary by a constant mean 
pressure gradient applied in the streamwise direction. Boundary conditions 
in the homogeneous directions are periodic whereas at the walls, no-slip and 
no-penetration conditions are imposed to the velocity field, and 9 is set to a 
constant value. A mean scalar gradient is thus established between the two 
walls, ensuring a mean passive scalar profile in the wall-normal direction. This 
bo undary condit i on for t he passive scalar has been alread y emp loyed for example 
by_ Kim &; Moinl (J1989I ). Ijohansson fc Wikstroml IJ1999I ). and iKawamura et al. 



( 20001 ). It is worth mentioning here that an alternative boundary condition can 



be used, that consists in imposing a constant scalar flu x at each wal l. Suc h 
a boundary condition has been chosen, for example, by Kasagi et al. 



iKasagi &: Ohtsubol (J1993I ) , IKawamura et all (|1998l ) and lKawamura et all ( 



1992), 



1999). 



The differences between the two boundary conditions have been addressed by 
Kawamura et all ( 20001) : they show that scalar statistics, such as the mean 



profile, the root-mean-square fluctuation and the turbulent flux, are basically 



unchanged in the near-wall region, while they differ in the core region due to 
the either zero or non-zero mean scalar gradient in the channel's mid plane. 

2.1 Validation 

The numerical code is validated by repli cating the DNS of a turbulent channel 
flow with passive scalar carried out by Uohansson fc Wikstroml ( 1999I ). This 



simulation is at i?e* = 265 and Sc = 0.71, with L x = 4tt5 and L z — 5.55, as in 
the original reference. The spatial resolution employed is N x = 256, N y = 256 
and N z = f92. After reaching a statistically steady state, our numerical simu- 
lation is continued for an overall averaging time of 4000 viscous time units, and 
60 statistically independent flow fields are writt en to disk for further an a lysis. 
Figure Q] compares the results obtained by Uohansson fe; Wikstroml (J1999J) 



to the present results. The two simulations are compared in terms of wall- 
normal profiles of the mean scalar (9) + , the root-mean-square value a^ of its 
fluctuations, and the mean scalar dissipation (eg) + . The superscript + indicates 
a non-dimensional quantity after scaling with wall (inner) variables. The inner 
scaling quantities for the relevant physical variables involved in this work are as 
follows. The inner velocity scale is given by the friction velocity u* = (j w / p) 1 / 2 , 
where t w and p are respectively the wall-shear stress and the fluid density, and 
the inner time scale £* is defined by £* = v/u 2 . For the passive scalar field, the 
inner scale 9* is given by: 

a _ 7 d(0> 



m* dy 

where the subscript w indicates a derivative evaluated at the wall. 

Given the definition ^ of the dissipation of scalar variance, its relevant 
scaling quantity e* is: 

7 fd{9) 



Sc \ dy 

A close inspection of figure[T]establishes the agreeme nt between the output of 



the pr esent numerical tool and the results obtained by Ijohansson fc Wikstrom 



(1999). The collapse of the curves is such that only for the profile of a^ can the 



two simulations be discerned. This little residual difference is most probably 
due to the different time span used for computing statistics. 

3 Computational parameters 

The present work focuses on the smallest scales of motion, and thus the available 
computational resources are best spent on spatial resolution. As a consequence, 
our simulations carried out to study the passive scalar field in turbulent plane 
channel flow have neither particularly high values of -Re* and Sc, nor particularly 
large sizes of the computational box. Its dimensions are set at L y = 26, L x — 
4.19J and L z = 2.095. The value of the friction Reynolds number is set at 
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Figure 1: Comparison between the present simulation (symbols) and results by 
Johansson fc Wikstroml (J19991) (dashed line) for i?e* = 265 and Sc = 0.71. All 
quantities in inner units. Left: wall-normal distribution of the mean passive 
scalar. Center: wall-normal profile of the r.m.s. value of the passive scalar fluc- 
tuations. Right: wall-normal mean profile of the dissipation of scalar variance. 
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Table 1: DNS parameters: grid spacings Aa;, Ay, Az are expressed in terms of 
(vb)w, the mean value of the Batchelor scale at the wall. In terms of viscous 
wall units, at the wall (t]b)w = 1-6. 



-Re* = 160. The value of the Schmidt number is set as 5c = 1, so that, according 
to Eq. ([1]), the mean Kolmogorov and Batchelor lengthscales are identical. 

Three DNS at increasing spatial resolution are carried out. They are labelled 
with L, M and H throughout the paper, to indicate Low, Medium and High 
resolution. The parameters defining the discretization of the three simulations 
are summarized in table Q] The grid spacings are expressed, as usual in wall 
turbulence, in viscous wall units, but they are also reported in terms of (t?b) w , 
the mean Batchelor length scale evaluated at the wall. 

The temporal discretization is adapted to the spatial discretization, so that 
the time step used in simulation H is the finest. The overall averaging time (2400 
viscous time units) and the number of statistically independent flow fields (60) 
stored on disk for further analysis for each of the three simulations are, however, 
left unchanged for all the simulations. 

The spatial resolution of simulation L is comparable to the resolutions em- 
ployed in most wall-turbulence DNS performed to date wit h passive scalars 
( Kawamura et all Il998l Il999t Ijohansson fc Wikstroml . [l999j). Such a spatial 



resolution is usually considered t o be adequate as far as the velocity field is 



concerned (jMoin fc Maheshl . 119981 ). 

On the other hand, simulation H has the highest resolution: owing to the 
increase with y of the lengthscale (77), that equals (t]b} since Sc = 1, the grid 
spacings Ax, Az and Ay are consistently smaller than (t/b) in the whole wall- 
normal range. Still, motions at scales smaller than (775) do exist. Figure [2] 
reports the statistical distribution of t/b measured at a few wall-normal loca- 
tions, and compares it to the spatial resolution of simulation H at the same y 
positions. It turns out that even the simulation employing the finest resolution 
could be only marginally resolved. At each wall-normal distance shown, in fact, 
either the streamwise (or spanwise) grid spacing, which is set constant to one 
wall unit throughout the channel, or the wall-normal spacing, that smoothly in- 
creases from the wall to the centerline, is larger than the smallest values locally 
assumed by 77^- 

The spatial resolution of simulation M is midway between L and H. Note 
that the streamwise and spanwise spatial resolution of case M, when expressed 
in terms of (j]b)w, is finer than the highest resolution employed in past DNS of 
passive scalar turbulent channel flows with Sc > 1 we are aware of. In absolute 



terms, this is the study recently carried out bv ISchwertfirm fe Manhartl ([20071 ) 
whose finest grid spacings are Ax + = 0.68, Az + = 0.85 and 0.18 < A?/ + < 0.75. 
Their DNS, however, was aimed at investigating flows with high values of the 
Schmidt number, so that the resulting spatial resolution evaluated in terms of 
the Batchelor scale is rather coarse: Ax/(t]b)w = 3.09, Az/(i]b) w — 3.86 and 
0.82 < Ay/{r) B ) w < 3.41. 

4 Statistical moments of eg 

The lowest-order statistical moments of the scalar dissipation eg are considered 
first. Figure [3] shows the wall- normal distribution of statistical moments of eg, 
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Figure 2: Probability density functions of 77^ for: (a) y + = 0; (b) y + = 5; (c) 
y + = 73; (d) y + = 160. The vertical lines indicate the spatial resolution of 
simulation H: the continuous vertical line is Ax + = Az + = 1, whereas Aj/ + is 
shown as a dashed vertical line, and increases with y. 
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Figure 3: Effects of the spatial resolution on the mean wall-normal distribution 
of the scalar dissipation rate eg (left) and on its variance a\ (right). 
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Figure 4: Effects of the spatial resolution on the mean wall-normal distribution 
of the turbulent kinetic energy dissipation rate e (left) and on its variance of 

(right). 
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respectively of order one (the mean profile) and order two (the variance), com- 
puted at different spatial resolutions. It can be easily appreciated how neither 
simulation L nor simulation M have the spatial resolution required to compute 
these integral quantities in a mesh-independent way. Indeed, a complete overlap 
between curve H and curve M cannot be observed, and thus even simulation 
H might still be slightly under resolved to correctly capture the full details of 
the fluctuations of eg . 

The three curves for (eg) + in figure [3] (left) collapse in the central region of 
the channel and in the logarithmic layer, indicating that in these regions case L 
has enough resolution to represent (eg) correctly. The three profiles are rather 
similar in the very-near-wall region too. Differences among the three curves 
become however evident in the buffer layer (say for 10 < y + < 30), where the 
separation becomes significant and a change of about 10% in the prediction of 
the mean value is observed. In this region, simulation H predicts lower absolute 
values of (eg), suggesting that on average the (absolute value of) dissipation of 
scalar variance associated to the finest scalar structures is lower. 

Spatial resolution effects appear to be more relevant when the variance of 
eg is considered: in figure [3] (right) the three curves remain separate for the 
whole y + range. The curves corresponding to cases L and M present a similar 
behavior in the viscous sub-layer, suggesting that, at least in this flow region and 
at these resolution levels, the wall-normal spatial resolution, which is identical 
in simulations L and M, could be the predominant factor. 

It is important to note how the same analysis carried out for the dissipation 
rate e of turbulent kinetic energy shows that these effects are peculiar to the eg 
field. Figure [3] should be confronted with figure 01 that shows the wall-normal 
profiles of the first two statistical moments of e. Though slight resolution effects 
can still be observed, particularly between cases L and M, it clearly emerges 
that the passive scalar field is characterized by a mu c h fine r and intermittent 
structure than the underlying velocity field ( Warhafti . 120001 ) . Hence turbulent 



flows with passive scalar, even at Sc = 1 and (tjb) = (??), present resolution 
requirements for correctly investigating integral-scale quantities of the scalar 
field which are more demanding than those required for the velocity field. When 
these requirements are not fully met, the insufficient resolution causes sizeable 
errors in the prediction of scalar-related mean quantities. 

5 PDF of e e 

In this Section, the resolution effects already observed in terms of integral quan- 
tities are discussed in terms of their intensity distribution. We consider first in 
i 35.ll particularly intense events, i.e. events characterized by extreme values of 
eg, and examine how the spatial resolution affects the right tail of the proba- 
bility density function of the quantity (eg ~~ (eg))/a eg . Then in §5.21 the same 
analysis is carried out for weak events, (i.e. events for which eg/ (eg) <C 1) by 
looking at the left tail of the probability density function of eg/ (eg). 
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5.1 Strong dissipative events 

Figure [5] shows the probability density function of the quantity z = (eg — 
{€g))/a eg computed at y + = (i.e. at the wall) and y + = 5 (in the viscous 
sublayer). Both the low- resolution simulations L and M are characterized by 
PDF with wide tails, which appear to shrink when the spatial resolution is in- 
creased in case H. This implies that, in the near-wall region, extremely intense 
events, that are present in cases L and M disappear in simulation H. This fea- 
ture becomes more evident in the viscous sublayer, where the largest difference 
among the tails of the three curves is observed. The limited spatial resolution of 
case L also results in overestimating events of intermediate intensity. Overall, 
these dissimilarities account for the higher absolute value of (eg) in the near-wall 
region previously observed in figure [3] for cases L and M. 

These conclusions, based on statistical quantities, can be arrived at by also 
observing instantaneous snapshots of the eg field. In particular, a wall-parallel 
section of the computational domain taken at y — (i.e. at the wall) is shown 
in figure [51 where the morphology of an eg field computed with the lowest (top) 
and the highest (bottom) resolutions are compared. The qualitative differences 
between the two slices help explaining the separation of the curves observed 
above in the context of figure [5j Case L always shows a larger number of regions 
with intense eg; their size is larger when compared to the analogous structures 
observed in case H. Furthermore, the typical values of eg for events of extreme 
intensity appears to be larger in simulation L. In this pair of snapshots, for 
example, the instantaneous maximum of e~g is 4.15 for case L and 2.98 for case 
H. Analogous characteristics can be observed in wall-parallel cuts taken through 
the viscous sublayer. 

The pattern just described, i.e. one where marginal resolution tends to 
overemphasize extreme dissipations, reverses when the attention is shifted from 
the near- wall region to the logarithmic layer and the central region of the chan- 
nel. Figure [7] reports again the PDF of eg, this time computed at y + = 73 
and y + — 160 (channel centerline). Here the most intense scalar dissipation 
events are reliably predicted by the highest spatial resolution only. This fea- 
ture is more evident in the log-layer, where the tails of the probability density 
functions corresponding to cases M and H are definitely longer than the tail of 
the curve for case L. Moreover, in the log-layer, the intensity of the strongest 
events captured by the simulation increases with the adopted spatial resolution, 
as evident from the length of the PDF tails. This implies that resolution M is 
certainly inadequate to capture the most intense scalar dissipation events in the 
log layer, and the same might be true for case H too. At the channel centerline, 
the separation between the PDF tails becomes smaller, yet case L is evidently 
not capable of capturing the strongest events. 

Again, a look at instantaneous fields of eg in wall-parallel planes confirms 
these remarks. Figure [5] shows sections at y + — 160 for simulations L and H, 
and shows how the resolution of case L does not represent the eg field correctly, 
owing to its extremely intermittent spatial structure. The maxima of e@ in these 
two particular slices are 1.65 (case L) and 4.35 (case H). 
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Figure 5: PDF of z = (e$ — (eg))/<T ee for two wall-normal positions at y + = 
(top) and y + = 5 (bottom). 
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Figure 6: Two-dimensional instantaneous wall-parallel cuts of the scalar dissi- 
pation field eg for y + = for simulation L (top) and H (bottom). Levels from 
by 0.25 increments. Maxima arc 4.15 (top) and 2.98 (bottom). 
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Figure 7: PDF of z = (eg — (eg))/a ee for two wall-normal positions at y + = 73 
(top) and y + — 160 (bottom). 
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Figure 8: Two-dimensional instantaneous wall-parallel cuts of the scalar dissi- 
pation field e^ for y + = 160 for simulation L (top) and H (bottom). Levels 
from by 0.25 increments. Maxima are 1.65 (top) and 4.35 (bottom). 
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Figure 9: PDF of z = ee/(e e ) for y+ = 5 (left), and y+ = 73 (right). 



5.2 Weak dissipative events 

Figure [9] shows the PDFs of the quantity z — eg /(eg) for the three different 
simulations in the near-wall region and in the log-layer. The choice of this 
quantity, and the logarithmic scale used on the horizontal axis, emphasize the 
left side of the PDF, and highlight the effects of the spatial resolution on the 
weakest dissipation events. 

The emerging pattern is that the left tail of the PDF for case L is con- 
sistently shorter than the tails for the other simulations. This means that a 
marginal spatial resolution is not capable to identify extremely we ak events of 
scalar dissipation. As already observed bv lSchumacher et all (J2005I ). this can be 
explained by the fact that numerical simulations with poor resolution are char- 
acterized by relatively high noise (when compared to more resolved simulations) 
and this implies a lower signal-to-noise ratio. The noise becomes particularly 
significant in the regions of low-magnitude eg events, that end up being covered 
by the noise floor. 

Simulations M and H, on the other hand, show no dissimilarities as far as 
the weakest eg events are considered, if exception is made for the viscous sub- 
layer (see again figure [9|), where the highest resolution is needed to observe the 
smallest events. 
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Figure 10: Three dimensional visualization of iso-scalar dissipation surfaces 
corresponding to eg/e* = 0.73 in the bottom half of the channel. 



Looking at the same quantities in the central region of the channel (not 
shown) confirms these observations, whereas at the wall, on the opposite, the 
three curves essentially coincide at small z. This last observation is explained 
both by the isoscalar boundary condition employed in the present work (which 
implies a tendential uniformity of the scalar in wall-parallel planes for very small 
wall distances) and by the larger characteristic length scale of the structures 
in the near-wall region, which can therefore be represented with lower spatial 
resolution. 

6 Discussion and conclusions 

In the present work three Direct Numerical Simulations of a passive scalar in 
turbulent channel flow have been carried out at different spatial resolutions. 
The aim was to investigate the effects of the spatial resolution on the statistical 
features of the flow, and to describe their change with the distance from the 
wall. We have centered our analysis on the statistical description of the scalar 
dissipation eg . Grid-dependent results have been observed not only in the statis- 
tical description of the passive scalar's micro-structure, but - more surprisingly 
- also in its integral-scale characteristics, unless the grid spacing in all the three 
directions is smaller than the local value of the Batchelor length scale. 
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The wall- normal profiles of (eg) and its variance show that, due to the strong- 
coupling between large-scale advective motions and small-scale molecular diffu- 
sive mixing, the most resolved DNS is needed even for the correct description 
of such integral-scale quantities. The required resolution goes far beyond what 
has been routinely adopted in the past for similar simulations. The resolution 
requirements become of course more stringent as soon as the order of the con- 
sidered statistical moments increases. Even though such a specific study has 
probably never been carried out for the velocity field, a comparison between the 
wall-normal mean profiles of eg and its variance with the corresponding profiles 
for the dissipation rate of turbulent kinetic energy strongly suggests that such 
strict requirements concern primarily the passive scalar, owing to its distinctive 
intermittent character. 

Extra-fine spatial resolution becomes essential when the smallest scales of 
the passive scalar are the main feature one is interested in. The PDF of eg 
fluctuations, measured in the channel's log-layer and the mid plane and reported 
in ij5.fl and £)5.2[ do in fact show that the range of scales which are not resolved 
by standard spatial resolutions is responsible for the extremely intermittent 
nature of the passive scalar field, that is observed for both very weak and very 
int ense scalar dissipa t ion ev ents. This is coherent wi t h the pattern described 



by ISchumacher et all (|2005h and IWatanabe fc Gotoh] (|2007|) for homogeneous 



isotropic turbulence. When the wall is approached, on the other hand, the 
influence of the solid boundary reverses the picture. In this region of the flow 
we have shown how the incorrect representation of the smallest scales leads to 
a significant overestimate of the strongest eg events. 

The varying effect of insufficient spatial resolution in the different layers 
of the channel flow is related to the changes of the eg field throughout the 
channel. The relatively large, high-dissipation structures which can be ob- 
served in near-wall horizontal sections of the flow field, visualized in figure [6j 
do disappear, in fact, as the distance from the wall increases, and are replaced 
by thin and elongated structures, shown in figu re [HI which recall, as it could 



have been expected, the structures observed by iBrethouwer et all ([20031 ) and 



ISchumacher et all ( 20051 ) in homogeneous and isotropic turbulence. According 
to these papers, the thin and elongated structures are twisted, folded, very close 
and almost parallel to each other, and are most likely cross-sections of the sheet- 
like structures of eg which can be observed in the three-dimensional visualization 
of the scalar dissipation field plotted in figure [10] The sheet-like structures are 
predominantly parallel to the wall in the near-wall region, and assume more 
of an isotropic character in the central region of the channel. The decreasing 
characteristic dimensions of high-eg structures for increasing distances from the 
wall explains why intense eg events are underestimated by simulation L in the 
channel's log-layer and near the mid plane of the channel. 

The (at least partial) grid-dependency presented by several results discussed 
throughout this paper implies that even the highest spatial resolution employed 
(simulation H) could be only marginally adequate. This is also suggested by 
the comparison, carried out in figure [2 between the statistical distribution of 
the Batchelor's scale ijb and the spatial resolution of case H in the different 
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layers of the channel flow. Grid resolutions finer than the smallest local value of 
rjs will therefore be required to investigate small-scale features of passive scalar 
mixing in wall-bounded turbulent flows. 
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